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Abstract. We have analysed high resolution adaptive optics (AO) science demonstration data of the young, massive 
stellar cluster Arches near the Galactic Center, obtained with the Gemini North telescope in combination with 
the University of Hawai'i AO system Hokupa'a. The AO H and K' photometry is calib rated using HST/NICMOS 



observations in the equivalent filters F160W and F205W obtained by Figer et al. ( 1999 ). The calibration procedure 
allows a detailed comparison of the ground-based adaptive optics observations against diffraction limited space- 
1 based photometry. The spatial resolution as well as the overall signal-to-noise ratio of the Gemini/Hokupa'a data 

is comparable to the HST/NICMOS data. The low Strehl ratio of only a few percent is the dominant limiting 
factor in the Gemini AO science demonstration data as opposed to space-based observations. After a thorough 
technical comparison, the Gemini and HST data are used in combination to study the spatial distribution of 
stellar masses in the Arches cluster. Arches is one of the densest young clusters known in the Milky Way, with 
a central density of ~ 3 • 10 s Mq pc -3 and a total mass of about 10 4 Mq. A strong colour gradient is observed 
|' over the cluster field. The visual extinction increases by AAv ~ 10 mag over a distance of 15" from the cluster 

1~ \ core. Extinction maps reveal a low-extinction cavity in the densest parts of Arches (R < 5"), indicating the 

depletion of dust due to stellar winds or photo-evaporation. We correct for the change in extinction over the field 
and show that the slope of the mass function is strongly influenced by the effects of differential extinction. We 
obtain present-day mass function slopes of F ~ —0.8 ± 0.2 in the mass range 6 < M < 65 Mq from both data 
sets. The spatial analysis reveals a steepening of the mass function slope from close to zero in the cluster center 
to r ~ —1.7 ± 0.7 at 7? > 10", in accordance with a Salpeter slope (r = —1.35). The bias in the mass function 
towards high-mass stars in the Arches center is a strong indication for mass segregation. The dynamical and 
relaxation timescales for Arches are estimated, and possible mass segregation effects are discussed with respect 
to cluster formation models. 

Key words. Open clusters and associations: individual: Arches - Stars: luminosity function, mass function - Stars: 
early-type - Stars: formation - ISM: dust, extinction - Instrumentation: adaptive optics 



1. Introduction 

The Galactic Center (GC) is the most extreme star form- 
ing environment within the Milky Way. High stellar and 
gas densities, turbulent motion, tidal torques exerted by 
the steep gravitational potential, magnetic fields and an 
intense radiation field determine the physical environment 
of star formation in the GC region. Although disruptive 
forces exerted by the gravitational and radiation fields 
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counteract the agglomeration of material, the high gas and 
dust densities cause star formation in the GC environment 
to be most efficient. In particular, the formation of high 
mass stars and massive clusters is more successful than in 
any other region of the Milky Way. 

A detailed study of star formation processes and the 
stellar content of the GC region has until recently been 
limited to the brightest and most massive stars due to 
the large amount of extinction (Ay ~ 30 mag) along 
the line of sight. Additional constraints are imposed due 
to the spatial resolution at the GC distance of ~ 8 kpc 
(DM = 14.47 ± 0.08 mag, e.g., McNamara et al. p000|), 
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much farther than nearby star forming regions such as 
the Orion or p Ophiuchi star forming complexes, which 
have been studied in greater detail to date. Only with the 
advent of deep, high resolution near-infrared instruments, 
the analysis of stellar populations in young star clusters 
near the GC has become feasible. 

During the past few years, it has become evident that 
three out of four young starburst clusters known in the 
Milky Way are located in the GC region - namely, the 
Arches and Quintuplet clusters, as well as the Galactic 
Center Cluster itself. With a cluster age of only a few 
Myr for Arches and Quintuplet, the question arises how 
many clusters do actually form in the densest environment 
of the Milky Way. The 2MASS database yielded new in- 
sights into the estimated number of star clusters hid den in 
the dense stellar background. Dutra & Bica ( 2000 2001 ) 
report the detection of new cluster candidates of various 
ages located in the innermost 200 pc of the Galaxy found 
in 2MASS. Numerical simulations by Portegies Zwart et 



al. ( 2001 ) suggest that clusters with properties similar to 
the massive Arches and Quintuplet may have formed in 
the past in the innermost 200 pc, but were then dispersed 
and are now indistinguishable from the dense stellar back- 
ground. As dynamical evolution timescales are short due 
to the strong tidal field in the GC region (Kim et al. 1999| ) , 
young star clusters are disrupted quickly after formation, 
contributing to the Galactic bulge population. Thus, only 
the youngest clusters remain intact for the study of star 
formation in this extraordinary environment. 

The Arches cluster, at a projected distance of only 25 
pc from the GC (assuming a heliocentric distance of 8 
kpc to the GC), is one of the most massive young clusters 
known in the Milky Way. With an estimated mass of about 
10 4 Mq and a central density of 3 • 10 5 Mq pc~ 3 , Arches 
is the densest young star cluster (YC) known (Figer et al. 
199S ). From physical properties of Wolf-Rayet stars, the 
age of the cluster is estimated to be between 2 and 4.5 



Myr (Blum et al. 2001 ) . The stellar content of Arches has 
been studied by Figer et al. ( |l999[) using HST/NICMOS 
data. They derived a shallow initial mass function in the 
range 6 < M < 120 Mq with a slope of T = -0.7±0.1, but 
with significant flattening observed in the innermost part 
of the cluster (r = -0.1 ± 0.2). 

Most young star clusters and associations in the Milky 
Way display a mass function close to a Salpeter (1955) 
power law with a slope of T = —1.35. Several such 
star forming regions have been studied by Massey et al. 
( 1995a ), yielding slopes in the range —0.7 < T < —1.7 
with an average of —1.1 ± 0.1, which leads these authors 
to conclude that within the statistical limits no deviation 
from a Salpeter slope is observed. 

A flat mass function as observed in Arches implies an 
overpopulation of the high-mass end as compared to "nor- 
mal" clusters. The special physical conditions in the GC 
region have been suggested to enhance the formation of 
massive stars, thereby resulting in a flattened mass func- 



lapse and subsequent accretion model, as radiation pres- 
sure from the growing star is capable of reversing the gas 
infall as soon as the mass is in excess of 10 Mq (Yorke & 
Kriigel 1977). Assuming disk accretion instead of spher- 
ical infall, the limiting mass may be increased to 15 Mq 
(Behrend & Maeder 2001), still far below the mass ob- 
served in O-type stars. Various scenarios are suggested to 
circumvent this problem. Simulations with enhanced ac- 
cretion rates and collision probabilities in dense cluster 



centers (Bonnell et al. 1998), as well as growing accretion 
rates depending on the mass of the accreting protostar 



(Behrend & Maeder |200l[ ), allow stars of up to 100 Mq to 
form in the densest regions of a rich star cluster. In case 
of the GC environment, a higher gas density may lead 
to a higher accretion rate and/or to a longer accretion 
process in the protostellar phase. As long as the gravi- 
tational potential is strongly influenced by the amount 
of gas associated with the cluster, gas infall causes a de- 
crease in cluster radius and subsequent increase in the 
collision rate, reinforcing the formation of high-mass stars. 
Physical processes such as gravitational collapse or cloud 
collisions scale with the square root of the local density, ^/p 
(Elmegreen |1999| |200lD , causing an enhanced star forma- 
tion rate (SFR) in high density environments. Elmegreen 
( |200l| ) shows that the total mass as well as the maximum 
stellar mass in a cluster strongly depends on the SFR and 
local density. This is confirmed by observations of high- 
mass stars found predominantly in the largest star forming 
clouds (Larson 1982). 

Both the growing accretion and the collision scenario 
predict the high-mass stars to form in the densest central 
region of a cluster, leading to primordial mass segregation, 
which may be evidenced in a flat mass function in the 
dense cluster center. As an additional physical constraint, 
both scenarios require the lower-mass stars to form first, 
and the highest-mass stars last in the cluster evolution 
process. As the strong UV-radiation field originating from 
hydrogen ignition in high-mass stars expells the remaining 
gas from the cluster center, the accretion process should 
be halted immediately after high-mass star formation. 

The short dynamical timescales of compact clusters 
are, however, influencing the spatial distribution of stel- 
lar masses as well. On the one hand, high-mass stars are 
dragged into the cluster center due to the gravitational po- 
tential of the young cluster. On the other hand, low-mass 
stars may easily be flung out of the cluster due to inter- 
action processes, especially given star densities as high as 
in the Arches cluster. The result of these processes would 
also be a flat mass function in the cluster center, steep- 
ening as one progresses outwards due to dynamical mass 
segregation. Dynamical segregation is predicted to occur 



tion (Morris 1993). The formation of high-mass stars in 



itself poses serious problems for the standard core col- 



within one relaxation time (Bonnell & Davies 1995), which 
for compact clusters is only one to a few Myr, and should 
thus be well observable in Arches in the form of a spatially 
varying mass function. 

In addition to the internal segregation process, the ex- 
ternal GC tidal field exerts shear forces tearing apart the 
cluster entity. N-body simulations by Kim et al. (2000) 
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yield tidal disruption timescales as short as 10 to 20 Myr 
in the GC tidal field. We expect to find a mixture of all 
these effects in the Arches cluster. 

We have analysed adaptive optics (AO) data ob- 
tained under excellent seeing conditions with the Gemini 
North 8m telescope in combination with the University 
of Hawai'i (UH) AO system Hokupa'a. We are investigat- 
ing the presence of radial variations in the mass distribu- 
tion within the Arches cluster. We compare our ground- 
based results in detail to t he HST/NICMOS data pre- 
sented in Figer et al. (199E, hereafter FKM), discussing 
possible achievements and limitations of ground-based, 
high-resolution adaptive optics versus space-based deep 
NIR photometry. 

In Section 2, we will introduce the data and describe 
the reduction and calibration processes. In this context, 
a thorough investigation of the quality of ground-based 
adaptive optics photometry as compared to space-based 
diffraction limited observations will be presented. The 
photometric results derived from colour-magnitude dia- 
grams and extinction maps will be discussed in Section 
3. A comparison of Gemini and HST luminosity functions 
will be given in Section 4. The mass functions will be 
derived in Section 5, and their spatial variation will be 
discussed with respect to cluster formation scenarios. We 
will estimate the relevant timescales for cluster evolution 
for the Arches cluster in Section 6, and discuss the im- 
plication on the dynamical evolution of Arches. We will 
summarise our results in Section 7. 



2. Observations and Data Reduction 

2.1. Gemini /Hokupa'a Data 

We analysed H and K' images of the Arches cluster center 
obtained in the course of the Gemini science demonstra- 
tion at the Gemini North 8m telescope located on Mauna 
Kea, Hawai'i, at an altitude of 4200 m above sea level. 
Gemini is an alt-azimuth-mounted telescope with a mono- 
lithic primary mirror and small secondary mirror opti- 
mised for IR observations. The telescope is always used 
in Cassegrain configuration with instruments occupying 
either the upward looking Cassegrain port or one of three 
sideward facing ports. The University of Hawai'i adaptive 
optics (AO) system Hokupa'a is a 35 e lement curvature 
sensing AO system (Graves et al. 2000 ), which typically 
delivers Strehl ratios between 5 % and 25 % in the K- 
band. 

Hokupa'a is operated with the near-infrared camera 



QUIRC (Hodapp et al. |1996| ), equipped with a 1024 x 1024 
pixel HgCdTe array. The plate scale is 19.98 milliarcsec- 
onds per pixel, yielding a field of view (FOV) of 20.2". The 
images are shown in Fig. 0, along with the HST/NICMOS 
images used for calibration. 

The observations were carried out between July 3 and 
30, 2000. 12 individual 60 s exposures, dithered in a 4 
position pattern with an offset of 16 pixels (0.32") be- 
tween subsequent frames, were coadded to an i7-band im- 



age with a total integration time of 720 s. In K', the set 
of 34 dithered 30 s exposures obtained under the best ob- 
serving conditions was coadded to yield a total exposure 
time of 1020 s. The full width at half maximum (FWHM) 
of the point spread function (PSF) was 9.5 pixels (0.19") 
in H and 6.8 pixels (0.135") in K' . The observations were 
carried out at an airmass of 1.5, the lowest airmass at 
which the Galactic Center can be observed from Hawai'i. 

The i?-band data were oversampled, and 2x2 binning 
was applied to improve the effective signal-to-noise ratio 
per resolution element, allowing a more precise PSF fit. 
A combination of long and short exposures has been used 
to increase the dynamic range. For the short exposures, 
3 frames with 1 s exposure time have been coadded 
in H, and 16 such frames in K'. See Table |l| for the 
observational details. In the long exposures, the limiting 
magnitudes were about 21 mag in H and 20 mag in 
K' . Note that the completeness limit in the crowded 
regions was significantly lower. The procedure used for 
comple teness correction will be described in detail in 
section 2.1.6. 



2.1.1. Data reduction 

The data reduction was carried out by the Gemini data 
reduction team, F. Rigaut, T. Davidge, R. Blum, and A. 
Cotera. The procedure as outlined in the science demon- 
stration report^] was as follows: Sky images, obtained after 
the short observation period when the Galactic Center was 
in culmination, were averaged using median clipping for 
star rejection, and then subtracted from the individual im- 
ages. The frames were then flatfielded and corrected for 
bad pixels and cosmic ray hits. After inspecting the in- 
dividual frames with respect to signal-to-noise ratio and 
resolution, and background adjustment, the images with 
sufficient quality were combined using sigma clipped aver- 
aging. The final images were scaled to counts per second. 
For the analysis presented in this paper, this set of images 
reduced by the Gemini reduction team has been used. 



2.1.2. Photometry 

The photometry was performed using the IRAF fl (Tody 



1993) DAOPHOT implementation (Stetson 1987). Due to 



the wavelength dependence of the adaptive optics correc- 



The description of the Galactic Center dataset ca n be 
found at bttp:/ /www. gemini.edu/gallcry /observing/ 



release_doc/manual.html, and the Gemini North science 
demonstration data are publicly a vailable 



at 



[http://www.gemmi.edu/sciops/data/dataSV.html . 

IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 



4 



A. Stolte et al.: Mass Function of the Arches Cluster 



Fig. 1 is available in jpeg format on astro-ph 



Fig. 1. Upper panels: Gemini/Hokupa'a H (left) and K' (right) images, 20"x 20"(0.8 pc x 0.8 pc) with a spatial 
resolution of less than 0. 2". N orth is up and East is to the right. Lower panels: HST/NICMOS F160W and F205W 
images from Figer et al. (1999), geometrically transformed to the Gemini FOV. 



Table 1. Gemini/Hokupa'a and HST/NICMOS observations 



Date 


Filter 


single exp. 


Hcxp 


exp. total 


det. limit 


Oback 


resolution 


diffraction limit 


Gemini 


05/07/2000 


H 


1 s 


3 


3 s 


18.5 mag 


7.74 


0.17" 


0.05" 


05/07/2000 


H 


60 s 


12 


720 s 


21 mag 


0.19 


0.20" 


0.05" 


30/07/2000 


K' 


1 s 


16 


16 s 


17.5 mag 


3.73 


0.12" 


0.07" 


09/07/2000 


K' 


30 s 


34 


1020 s 


20 mag 


0.22 


0.13" 


0.07" 


HST 


14/09/1997 


F160W 






256 s 


21 mag 


0.04 


0.18" 


0.17" 


14/09/1997 


F205W 






256 s 


20 mag 


0.15 


0.22" 


0.21" 



tion and anisoplanatism^j over the field, the H and K' 
data have been treated differently for PSF fitting. While 
in H the PSF radius increases significantly with distance 
from the guide star, with a radially varying FWHM in 
the range of 0.18" to 0.23" (see the science demonstration 
data description), the K 1 PSF was nearly constant over 
the field (0.125" to 0.135"). This behaviour is expected 
from an AO system, as the isoplanatic angleQ 9q varies 
as A 6 / 5 , yielding a 1.4 times larger #0 in K' than in H, 
resulting in a more uniform PSF in K' across the field of 
view. 

As obscuration due to extinction decreases with in- 
creasing wavelength, many more stars are detected in K' 
than in H . For comparison, the number of objects found 
with K' < 20 mag and uncertainty <jk' < 0.2 mag was 
1017 (1020 s effective exposure time), while for H < 20 
mag and an < 0.2 mag we detected only 391 objects (720 
s effective exposure time), where in both cases visual in- 
spection led to the conclusion that objects with photo- 
metric uncertainties below 0.2 mag were real detections. 
On the other hand, the increased stellar number density 
in K' leads to increased crowding effects, such that we 
decided to use a non-variable PSF for the if' -band data 
after thorough investigation of the results of a quadrati- 
cally, linearly or non-varying PSF. It turns out that, due 
to a lack of isolated stars for the determination of the 
PSF variation across the field, the mean uncertainty is 
lower and the number of outliers with unacceptably large 



3 The anisoplanatism describes the degradation of the AO 
correction due to differences between wave fronts coming from 
different directions. 

4 The isoplanatic angle is defined as the angular distance for 
which the rms error of the wave front is less than 1 radian. 



uncertainties reduced when a non-variable PSF is used. 
Thus, the 5 most isolated stars on the K' image, which 
were well spread out over the field, were used to derive 
the median averaged PSF of the long exposure. In the 
case of the short exposure, where due to the very short 
integration time faint stars are indistinguishable from the 
background, leaving more 'uncrowded' stars to derive the 
shape of the PSF, 7 isolated stars could be used. In K' , 
the best fitting function was an elliptical Moffat-function 
with = 2.5. 

Due to the lower detection rate on the -ff-band image 
crowding is less severe, while the PSF exhibits more pro- 
nounced spatial variations than in K' . We thus used the 
quadratically variable option of the DAOPHOT psf and 
allstar tasks for our iJ-band images, with 27 stars to de- 
termine a median averaged PSF function and residuals. 
The best fitting function was a Lorentz function on the 
binned ii-band image. In both filters, the average FWHM 
of the data has been used as the PSF fitting radius, i.e., 
the kernel of the best-fitting PSF function, to derive PSF 
magnitudes of the stars. 

The short exposures have been used to obtain pho- 
tometry of the brightest stars, which are saturated in 
the long exposures. The photometry of the long and 
short exposures agreed well after atmospheric extinction 
correction in the form of a constant offset had been 
applied. The saturation limit was 13.0 mag in H and 13.3 
mag in K' . At fainter magnitudes, the photometry of both 
images was indiscernible within the uncertainties for the 
bright stars, and the better quality long exposure values 
were used. Furthermore, a comparison of the bright 
star photometries was used to estimate photometric 
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Fig. 2. Formal DAOPHOT photometric uncertainties. Top panel: Gemini/Hokupa'a H and K 1 photometry trans- 
formed to HST/NICMOS F160W and F205W filter magnitudes, ml60 and m205. Bottom panel: HST/NICMOS 
m!60, m205 photometry of the same 20 x 20" field (shown in Fig. fy). 



uncertainties (see Sect. 2.1.6 and Table 



2.1.3. Photometric Calibration 

To transform instrumental into apparent magnitudes, we 



used the HST/NICMOS photometry of Figer et al. ( |1999[ ) 
as local standards. The advantage of this procedure lies 
in the possibility to correct for remaining PSF deviations 
over the field, e.g., due to a change in the Strehl ratio 
with distance from the guide star or due to the increased 
background from bright star halos in the cluster center. 
Indeed, as will be discussed below, the spatial distribution 
of photometric residuals shows a mixture of these effects. 

We were able to use approximately 380 stars to de- 
rive colour equations. The residuals obtained for these 
stars after calibration allow a detailed analysis and cor- 
rection of field variations. The colour equations to trans- 
form Gemini instrumental H and K' magnitudes to mag- 
nitudes in the HST/NICMOS filter system were deter- 
mined using the IRAF PHOTCAL package, yielding: 



ml60 = H h 



m205 = K{ nst 



0.001(±0.017)-(if-iOin S t 

-0.028(±0.031) mag 
■ 0.023(±0.008) • {H - K ') inst 

-1.481(±0.016) mag, 

where i?i ns t and K[ nst are the Gemini instrumental mag- 
nitudes, and ml60 and to205 correspond to magnitudes 
obtained with the NICMOS broadband filters F160W and 
F205W, respectively. After the transformation had been 



applied, it turned out that the residual magnitudes in the 
two independent fitting parameters H — K' and K' still 
varied systematically over the field. As can be seen in Fig. 
^, the variation is not a simple radial variation increasing 
with distance from the guide star (GS), but a mixture of 
displacement from the GS position (shown in the contour 
map in Fig. || as a cross), and the position with respect to 
the cluster center or bright stars in the field. Fortunately, 
the variation was well behaved in the Indirection, and 
could be fitted by a fourth order polynomial. Close inspec- 
tion showed that two areas on the QUIRC array showed a 
remaining photometric offset compared to the HST pho- 
tometry. In the region 400 < X < 700 pixels and Y < 250 
pixels, the magnitude was underestimated by 0.1 mag. For 
X > 850 pixels and Y > 850 pixels, i.e., the upper right 
corner, the K' magnitude was overestimated by 0.25 mag 
(however, note that there are only 7 stars i n thi s region). 
For a discussion of these effects, see Sect. 2.1.4. We cor- 



rected for these offsets before deriving the F-correction, 
which then showed remarkable homogeneity over the en- 
tire field. This smooth correction function is probably due 
to discrepancies in the dome flat field illumination versus 
sky exposures. The correction was then applied to trans- 
formed K' and H — K' magnitudes, denoted -K^ rans and 
(H — K' )trans from now on. The .fftrans magnitude was 
calculated from the corrected if t 'rans an d (H — K') tranB 
values. 

An additional advantage of this procedure is the 
independence on uncertainties in colour transformations 
at large reddening and non-main sequence colours, as 
opposed to colour transformations derived from typical 
main sequence standard stars. Using the HST photometry 
as local standards, we are naturally in equal colour and 
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Fig. 3. Map of residuals of NICMOS vs. Gemini photometry (orientation as in Fig. [I]). 
Left panels: m205 — K[ rans , right panels: ml60 — -fftrans 

Stars have been binned into areas of 50 x 50 pixels, and the value displayed shows the average of all stars in each 
bin. Statistical fluctuations are large due to the varying number of stars in each pixel intervall, but the overall trends 
are clearly visible. The position of the guide star is marked by a cross, the cluster center as determined from the 
HST/NICMOS F205W image as an asterisk, and stars denote stars resolved in the Gemini images with magnitudes 
brighter than K' ir&ua = 13 mag. The strong correlation between the position of bright stars and positive residuals 
reveals the tendency of PSF fitting photometry applied to AO data to overestimate the flux of bright sources, and 
underestimate the flux of faint sources in their vicinity (see text for discussion). 



temperature regimes, allowing the direct comparison 
of the Gemini and HST photometry. For most parts 
of the paper, we remain in the HST/NICMOS system. 
We use t he colour equations obtained in Brandner et 
al. ( 2001 , hereafter BGB) to transform typical main 



diagram (Sect. |3.4[) and the derivation of the extinction 



sequence colours and theoretical isochrone magnitudes 
into the HST/NICMOS system where indicated. This 
allows us to transform mainly unreddened main-sequence 
stars, for which the BGB colour transformations have 
been established. The only exceptions are the two-colour 



variation from colour gradients (Sect. 3.1), where the 
extinction law is needed to determine the reddening 
path. We will use the notation "ml60" and "m205" as in 
FKM for magnitudes in the HST/NICMOS filters, and 
"fftrans" and "-ft^rans" f° r the Gemini/Hokupa'a data 
calibrated to the NICMOS photometric system. HST 
magnitudes transformed to the ground-based 2MASS 
system will be denoted by ( JH K s )2mass or simply JHK S . 
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2.1.4. Discussion of the residual maps 

The behaviour of the residual of the HST/NICMOS vs. 
Gemini magnitudes, mNicMOS-^Gcmini.trans, can be anal- 
ysed in more detail when studying the residual map and 
the smoothed contour plot. In (Fig. ||), positive (negative) 
residuals correspond to overestimated (underestimated) 
flux. From the map we denote a general tendency to over- 
estimate the flux. The contour maps show that positive 
residuals are correlated with the position of bright stars 
on the K' image, both in the crowded cluster center as well 
as in the area to the lower left, where a band of bright stars 
is located (see Fig. |). This is the area where the K' mag 



nitudes were found to be underestimated in Sect. 2.1.3 , 
and thus the flux overestimated. This suggests that the 
increased background due to the uncompensated seeing 
halos of bright stars (cf. Sect. 2.1.5), inherent to AO ob- 



servations, causes an overcstimation of the flux of bright 
(K' ~ 13 mag)) sources. On the other hand, points with 
negative residuals are mainly correlated with fainter stars 
(K ~ 16 mag), suggesting that this enhanced background 
leads to an oversubtraction of the individually calculated 
background of nearby fainter stars. The result is an un- 
derestimate of the flux of companion stars in the vicin- 
ity of bright stars. The correlation of positive residuals 
with the position of bright stars seems to be less pro- 
nounced in the -ff-band image (Fig. ||). In H, the distance 
from the guide star is supposed to be more important 
due to the smaller size of the isoplanatic angle at shorter 
wavelengths and consequently more pronounced anisopla- 
natism. Indeed, the smoothed residual contour plot shows 
a symmetry in the residuals around the guide star, with 
close-to-zero residuals in the immediate vicinity of the 
guide star, where the best adaptive optics correction can 
be achieved. With increasing distance from the guide star, 
the residuals increase not only towards the bright cluster, 
but also to the west (left in Fig. 1) of the field, indicating 
that remaining distortion effects from the AO correction 
are mixed with the problem of the proximity to bright 
stars as seen in K' . 

2.1.5. Strehl ratio 

The Strehl ratio (SR) is defined to be the ratio of the 
observed peak-to-total flux ratio to the peak-to-total flux 
ratio of a perfect diffraction limited optical system. This 
definition allows to compare the quality and photometric 
resolution of different optical systems using a single char- 
acteristic quantity. 

SR = (i ? peak/-Ptotal)obs/(-Ppeak/-Ptotal)theo 

where -F pC ak is the maximum flux value of the PSF, and 
-Ftotai is the total flux, including the uncompensated halos 
induced by the natural seeing. The labels obs and theo re- 
fer to the observed and the theoretically expected flux, re- 
spectively. Diffraction limited theoretical PSFs have been 
calculated using the imgen task of the ESO data analy- 



the Gemini ff-band 720s exposure is found to be 2.5 % 
compared to 95 % in the HST F160W image, and 7 % 
in K' (1020s) compared to - 90 % in F205W. The low 
SRs measured in the Gemini science demonstration data 
indicate that more than 90 % of the light of a star is dis- 
tributed into the resolution pattern of the AO PSF and 
the halo around each star induced by the natural seeing. 
It has turned out that these halos cause a significant lim- 
itation to the resolution and depth of the observations in 
a crowded field, as faint stars can be lost to the enhanced 
background in the vicinity of bright objects. When com- 
paring the HST and Gemini luminosity functions (Sect. ^), 
this effect causes the main difference between both data 
sets. 

In the case of very low Strehl ratios, the SR does not 
directly indicate the fraction of the flux concentrated in 
the FWHM area of the PSF. A much larger fraction of 
the source flux can be used for PSF fitting in this case, 
although the spatial resolution is limited by the large 
FWHM as compared to diffraction limited observations 
(cf. Tab. |). The ratio of the flux in the FWHM kernel of 
the compensated stellar image to the total flux, 



FR h 



"FWHM/^total 



may be determined by creating curves of growth for indi- 
vidual stars (Stetson 1990 ). The larger number of nearly 
isolated stars found on the i7-band image (used for PSF 
creation) allowed a reliable determination of Fi? Q b s only 
in H , although many stars on the If-band image were still 
too influenced by neighbours to study the aperture curve 
of growth in detail. We were able to create well-behaved 
curves of growth for 7 stars. As in the PSF fitting routine, 
we have used the FWHM of the PSF as kernel radius 
and as the reference flux for the flux ratio determination. 
These ratios range from 0.47 to 0.58, with an average of 
0.53 ± 0.05, i.e., ~ 50 % of the integrated point source 
flux are used for PSF fitting. In addition, the variation of 
the flux ratio Fi? Q b s over the field serves as an indicator of 
the Strehl ratio variation. If the AO correction mechanism 
is the dominant factor determining the concentration of 
the flux into the FWHM kernel, the Strehl ratio and thus 
Fi? bs should decrease with distance form the guide star, 
as the seeing correction worsens. No correlation of the flux 
ratio with distance from the guide star is found. Though 
these are small number statistics, this supports the sug- 
gestion that the sensitivity variations over the field are not 
predominantly due to increasing distance from the guide 
star (Sect. |2.1.4[ ). 



sis package eclipse (Devillard 1997). The Strehl ratio in 



2.1.6. Photometric uncertainties and incompleteness 
calculation 

For the incompleteness correction, artificial frames 
were created with randomly positioned artificial stars. 
Magnitudes were also assigned automatically in a random 
way. Due to the very crowded field, only 40 stars were 
added to each artificial frame in order to avoid significant 
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Fig. 4. Incompleteness tests performed on the Gemini H 
and K' data. The plot demonstrates the necessity to in- 
clude the matching between single filter observations to 
obtain a realistic estimate on the incompleteness, when 
matching of objects is needed to create colour-magnitude- 
and two-colour-diagrams. 



changes in the stellar density. A total of 100 frames was 
created for both the H and K 1 deep exposures, leading to 
a total of 4000 artificial stars used in the statistics. In ad- 
dition to the individual incompleteness in each band, the 
loss of sources due to scatter of the main sequence gener- 
ated by the more uncertain photometry in the dense parts 
of the cluster was estimated. For this purpose the artifi- 
cial K' stars were assigned a formal instrumental 'colour' 
of (H — K') = 0.33 mag (corresponding to 1.745 mag 
after photometric transformation), derived from the av- 
erage observed instrumental colour of the main sequence, 
and via this transformed into instrumental H magnitudes. 
Artificial stars were inserted at the same positions in the 
H and K' frames. In this way the procedure also accounts 
for stars lost due to the matching of H and K 1 data. The 
artificial stars were calibrated using the colour equations 
shown in Sect. 2.1.3 , thus allowing us to estimate the loss 
of stars in the mass function derivation due to the applied 
main sequence colour selection (see Sect. ||). This resulted 
in significantly larger corrections as compared to the in- 
dividual filter recoverage without matching and colour se- 
lection. As an example, the results for the mass function 
calculation performed on the artificial stars are displayed 
in Fig. |. 

As the recovery rate depends strongly on the stellar 
density and thus radial distance from the cluster center, 
the incompleteness correction will be determined in depen- 
dence of the radial bin analysed when radial variations in 



the MF are studied (Sect. |5.3| ). 

In addition to luminosity and mass function correc- 
tions, the artificial star tests were used to estimate the real 
photometric uncertainties by comparing inserted to recov- 
ered magnitudes of the artificial stars. The median differ- 
ence between the original and the recovered magnitude 
of the artificial stars, Am art i = <m addc d - wdaophot>, 



has been used as an estimate of the real photometric un- 
certainty. To obtain the median uncertainty, the intcrvalls 
< Am art i < 1 and — 1 < Am art i < have been treated 
individually, and the mean of the absolute value of both 
median values, weighted with the number of objects in 
each intervall, is quoted in Tab. [| The overall flux devi- 
ation, including positive and negative deviations, is close 
to zero for stars brighter than 20 mag (< ±0.004), and for 
fainter stars becomes —0.17 and —0.13 in H and K' , re- 
spectively, showing a tendency to underestimate the flux 
of faint stars. This tendency is more serious in the cluster 
center, where comparably large uncertainties are already 
observed for magnitudes fainter than 18 mag in both path- 
bands. 

In a second test, photometry of the bright stars in the 
short exposures was compared to the magnitudes of the 
deep exposures, yielding Am s i = <m s hort — wiong> using 
the same procedure as for artificial stars. Note, however, 
that the quality of the short exposures is much worse than 
the one of the long exposures due to the high background 
noise. Therefore, the artificial star experiments, for which 
only the deep exposures have been used, yield a more real- 
istic estimate of the photometric uncertainties. The results 
of both tests are summarised in Table |[ 

The resulting uncertainty is roughly a factor of 2 to 3 
larger than the theoretical magnitude uncertainty ctdao 
determined from detector characteristics by DAOPHOT. 
The photometry in the cluster center shows a larger uncer- 
tainty than the photometry in the outskirts. As expected 
in a crowding-limited field, this also implies a reduced de- 
tection probability of faint sources in the center of the 
cluster. 

2.2. HST/NICMOS data 

HST/NICMOS observations have been obtained in the 
three broad-band filters F110W, F160W and F205W, 
roughly equivalent to J, H and K. The basic parame- 
ters are included in Table [l]. For a detailed description of 
the HST data and their reduction see Figer et al. (1999). 



3. Photometric results 

3.1. Radial colour gradient 

During the process of calibration we realised that a strong 
colour gradient is present in the Arches field. In Fig. 0, 
the H — K' colour is plotted against radial distance from 
the cluster center. As the extinction law has not yet been 
derived for HST filters, we have used the colour transfor- 
mations in BGB to transform NICMOS into 2MASS mag- 
nitudes. Though the 2MASS filters deviate slightly from 
the standard Johnson JHK filters used to determine the 



extinction law (Rieke & Lebofsky 1985), we will be able to 
estimate the approximate amount of change in visual ex- 
tinction across the field. The extinction parameters from 
Rieke & Lebofsky (1985) are given by 



Aj/A v = 0.282 A H /A V = 0.175 A K /A V = 0.112, 
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Table 2. Photometric uncertainties derived from artificial star experiments (Am art i) and the comparison of short and 
long exposures (Am s i). The photometric uncertainties determined by DAOPHOT are given for comparison. Note that 
for the brightest bin, 10 - 12 mag, only the photometry of the short exposures was available. The higher starting bin 
in H is due to the calibration procedure of the inserted artificial stars. Magnitudes are given in the NICMOS system 
(F160W, F205W). 



Band 


mag 


Am ar ti 


Am s i 


fDAO 


Am al . ti 


Am s i 


0"DAO 




Am s i 


0"DAO 






all 




R < 10" 


R > 10" 








H 


12 - 14 


0.007 


0.059 


0.004 


0.012 


0.052 


0.005 


0.004 


- 


0.002 




14 - 16 


0.015 


0.061 


0.008 


0.025 


0.065 


0.010 


0.011 


0.035 


0.004 




16 - 18 


0.044 


0.128 


0.015 


0.077 


0.120 


0.017 


0.038 


0.134 


0.009 




18 - 20 


0.119 




0.036 


0.167 




0.039 


0.098 




0.030 




> 20 


0.264 




0.142 


0.514 




0.142 


0.265 




0.144 


K' 


10 - 12 


0.004 




0.005 


0.004 




0.005 


0.003 




0.006 




12 - 14 


0.008 


0.042 


0.005 


0.011 


0.048 


0.005 


0.007 


0.028 


0.005 




14 - 16 


0.027 


0.072 


0.007 


0.041 


0.086 


0.007 


0.020 


0.042 


0.006 




16 - 18 


0.071 


0.121 


0.016 


0.125 


0.166 


0.017 


0.057 


0.088 


0.016 




18 - 20 


0.206 




0.060 


0.280 




0.073 


0.189 




0.056 




> 20 


0.411 




0.274 


0.520 






0.360 




0.274 



where Aputer 1S the extinction in the given filter. For a 
change in extinction of AAy this leads to 

AA V = A(Aj - A H )/0.W7 
AAy = A(A H - A K )/0.063 
AA V — A(Aj — Ak-)/0.170 

We are able to measure the right-hand side of each expres- 
sion from the transformed HST/NICMOS JHK S photom- 
etry. The resulting AAy is given in each plot in Fig. ^|, 
together with the fitting uncertainty from the rms scatter 
in the colour. 

From the linear fits in Fig. ^| we see that Ay increases 
by about one order of magnitude over the entire field 
when moving outwards from the cluster center. The ef- 
fect is most pronounced in the HST J — K s vs. radius 
diagram (Fig. ||[ bottom) , where the longest colour base- 
line is used. We derive a change in visual extinction of 
AAy = 10.71 ± 2.47 mag over the Gemini field (1000 pix- 
els = 0.8 pc). Notably, if only the innermost 5" (250 pixels, 
0.2 pc) are fitted, no variation in Ay is observed. When 
fitting the core separately, we get AAy — 0.77 ± 1.12 mag 
for R < 5" (250 pixels) versus AAy = 7.87 ± 2.85 mag for 
5" < R < 20" (250 < R < 1000 pixels). The latter value 
corresponds to AAy = 10.5 mag/ 1000 pix~ , consistent 
with the trend over the entire field. The small radial trend 
and low extinction value in the cluster center indicates the 
local depletion of dust. This could be either due to winds 
from massive stars or due to photo-evaporation of dust 
grains caused by the intense UV-radiation field. 

A change in Ay of ~ 10 mag/1000 pix" 1 is also 
consistent with the result found in J — H, 8.8 ± 
2.1 mag/ 1000 pix -1 , while a larger value of AAy = 
14.9 ± 3.2 mag is derived from the H — K s plot. Due to 
the uncertainties we conclude that the extinction varies 
by AAy ^9 — 15 mag across the Arches field of 20" or 
0.8 pc, most likely closer to the lower value. This is in 



any case a tremendous change in the dust column density 
along the line of sight, with strong implications on limiting 
magnitudes and the potential detection of faint objects. 

We have used a linear fit to the colour variation with 
radius for R > 5" to correct for the strong change in red- 
dening observed in the outer cluster field. The values for 
cluster center stars (R < 5") have been left unchanged, 
due to the large scatter and the very small trend found. 
Thus, these adjusted colours are scaled to the cluster cen- 
ter, where Ay is lowest. From the Rieke & Lebofsky (1985) 
extinction law, the change in if-band magnitude with ra- 
dius corresponding to the change in colour can be derived 
as AA K = 0.112/0.063 • AAh-k- We have used this re- 
lation to adjust the .fT-band magnitudes accordingly. The 
'dereddened' colour-magnitude diagrams will be shown in 
direct comparison with the observed CMDs in Sect. 3.3 
(cf. Fig. |). 

3.2. Extinction maps 

We can calculate the individual extinction for each ob- 
ject by shifting it along the reddening vector onto the 
isochronc. If we assume that intrinsic reddening plays only 
a minor role for most of the stars within the cluster, we can 
combine the individual reddenings to estimate the spatial 
extinction variation. In practice, we have applied a redden- 
ing of Ay = 15 mag to a 2 Myr isochrone of the Geneva set 
of models (Lejeune & Schaerer [2001 ). This choice of red- 
dening ensures that the isochrone serves as a blue envelope 
for the bulk of the cluster stars, which are significantly 
more reddened. The choice of the isochrone will be dis- 
cussed in context with the mass function (Sect. |5|), where 
the physical effects of the population model used are more 
important. As discussed in Grebel et al. ( 1996] ), colours of 
main-sequence stars also depend on binarity, stellar rota- 
tion, and intrinsic infrared excess. We ignore these effects 
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Fig. 5. Colour variation across the Gemini field, 
as observed in instrumental Gemini magnitudes and 
HST/NICMOS data. As a trend of increasing colour ex- 
cess with increasing distance from the Arches cluster cen- 
ter is observed in the two independent datasets, an in- 
strumental effect as the cause for this variation is highly 
improbable. 



here since we have no means to distinguish them from 
reddening effects. In the worst case, we will overestimate 
the reddening for some stars. We have then shifted the 
stars along the reddening path in colour-magnitude space, 
again assuming that the relative slopes of the extinction 
law approximately follow a Rieke & Lebofsky ( 1985| ) law. 
The resulting extinction map for the Gemini photometry 
is shown in Fig. [?]. 

Since the original HST/NICMOS field has twice the 
area of the Gemini central Arches field, we have also cal- 
culated the extinction map for the entire NICMOS field 
following the same procedure. The corresponding map is 
shown in Fig. |^. 

The extinction measured with this procedure in the 
JT-band lies in the range 1.9 < Ak < 4.1 mag, with an 
average value of 3.1 mag, correspondi ng to 16 < Ay < 37 
mag, <A V >= 27.7 mag. Cotera et al. ( 2000 ) derive a near- 
infrared extinction of 2.8 < Ak < 4.2 mag, with an aver- 
age value of <Ak> = 3.3 mag for 15 lines of sight towards 
several Galactic Center regions, corresponding to an aver- 
age visual extinction of < Ay> = 29.5 mag (transformed 
using Rieke & Lebofsky ( |1985| )). They obtain the high- 
est extinction towards a field close to the Arches cluster, 
Ay — 37.5 mag. This is very close to our highest extinc- 
tion value. The average value determined from individual 
dereddening here is the same as the average extinction ob- 
tained by Figer et al. fll999| ), <A V > = 27.7 mag. Note that 
the typical random scatter <j(Ak) from foreground dust 
density fluctuations found in GC fields is linearly related 
to the average extinction within a field. The relation de- 



termined by Frogel et al. (1999) from giant branch stars in 
22 pointings towards fields within 4° from the GC is given 
by a(A K ) = 0.056(±0.005) <A K > + 0.043(±0.005). This 
yields an expected natural scatter from GC clouds of only 
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Fig. 6. Colour trends over the Arches field as observed in 
the HST/NICMOS data set within the area covered by the 
Gemini observations. The radial distance from the cluster 
center is given in pixels on the Gemini scale. For the cal- 
culation of A Ay, the Rieke & Lebofsky (1985) reddening 
law has been assumed. 



(j(Ak ) — 0.22 mag for <Ak> = 3.1 mag, much below the 
difference in reddening observed in the Arches field. Thus, 
the change in extinction cannot be explained by the natu- 
ral fluctuations of the dust distribution in the GC region. 

Comparison of the cluster center main sequence pop- 
ulation with the main sequence colour of a theoretical 2 
Myr isoc hrone from the Geneva set of models (Lejeune & 
Schaerer 2001 ) , later-on used for the derivation of the mass 
function, yields an average extinction of Ay = 24.1 ± 0.8 
mag in the cluster center. This extinction value has been 
used to transform isochrone magnitudes and colours into 
the cluster magnitude system. It has been suggested that 
the brightest and most massive stars in Arches are Wolf- 



Rayet stars of type WN7 (Cotera et al. 199£; Blum et al. 
[2001 ). Fundamental parameters of Wolf-Rayet stars are 
compiled in Crowther et al. ( 1995| ). For stars of subtype 
WN7 they find typical colours of (H — K) ~ 0.2 mag, 
leading to an extinction of Ay = 24.9 ± 2.4 mag with 
an observed H — K colour of ~ 1.77 mag for the WN7 
stars, which were identified by comparison with the Blum 
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et al. (2001) narrow band photometry. This value is in 
very good agreement with the Ay determined from the 
main sequence colour in the cluster center. 



3.3. Colour-magnitude diagrams 

The resulting colour-magnitude diagrams for Gemini and 
HST are presented in Fig. ^ (upper panel). Two impor- 
tant differences are seen when inspecting the CMDs. First, 
the scatter in the main sequence is significantly larger in 
the ground-based photometry. While the HST/NICMOS 
CMD reveals a narrow main sequence in the cluster cen- 
ter (circles in Fig. the same stars display a much 
larger colour range in the Gemini CMD. The poor Strehl 
ratio in the Gemini/Hokupa'a data as compared to the 
HST/NICMOS data (see Sect. pX5| ) causes a high, non- 
uniform additional background due to uncompensated see- 
ing halos around bright stars, which decreases photomet- 
ric accuracy. In the dense regions of the cluster center, 
where crowding problems are most severe, the photom- 
etry is most affected. The number of faint, unresolved 
companion stars that merge into the high stellar back- 
ground underneath the bright cluster population is very 
high. As discussed in Sect. 2.1.4 , the halos of the bright 
stars hinder the detection of faint objects despite the prin- 
cipally high spatial resolution seen in individual PSF ker- 
nels. Operating at the diffraction limit, NICMOS is not 
restricted by these effects, yielding a better effective reso- 
lution especially in the dense regions. A tighter main se- 
quence and less scatter is the consequence. In Fig. ||, the 
innermost 5" of the Arches cluster are marked by open 
circles. It is clearly seen that most massive (bright) stars 
are located in the cluster center. 



A second effect observed is the much larger number 
of faint objects seen in the HST data (cf. Fig. |). As the 
limiting magnitude and the measured spatial resolution of 
the images are similar in both datasets, this, too, has to 
be a consequence of the low Strehl ratio in the AO data. 

The lower panel of Fig. ^ shows the "dereddened" 
CMDs, corrected for the radial colour gradient found an d 
the corresponding change in extinction, AAk (Sect. 3.1). 
The colours of stars beyond R > 5" have been adjusted 
to the colour of the cluster center. Comparison with the 
original CMDs shows that most of the bright, seemingly 
reddened stars fall onto the same main sequence after cor- 
recting for the colour trend. These stars are located at 
larger distances from the cluster center and thus suffer 
from more reddening by residual dust. As will be discussed 
in the context of the mass function (Sect. ||), these stars 
might have formed close to the cluster at a similar time 
as the cluster population. At the faint end of the CMD, 
there are, however, a large number of objects that re- 
main unusually red after the correction has been applied. 
These objects may either be pre-main sequence stars or 
faint background sources. Unfortunately, we are not able 
to disentangle these two possible contributions, and will 
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Fig. 9. Colour-magnitude diagrams. 
Left: Gemini/Hokupa'a, right: HST/NICMOS 
Lower pa nel: CMDs corrected for radial reddening gradi- 
ent (Sect. |t|) 



thus exclude objects significantly reddened relative to the 
main sequence when deriving the mass function. 

3.4. HST/NICMOS colour-colour diagram 

For comparison with the reddening path and a main se- 
quence in standard colours, the NICMOS filters have been 
transformed into the 2MASS JHK S system. In Fig. [l0[ 
we show the transformed HST/NICMOS colour-colour 
diagram for the stars bright enough to be observed in 
all three filters. The Ay values are from the Rieke & 



Lebofsky (1985) extinction law for standard JHK pho- 
tometry. Though we are aware of the uncertainties in- 
herent to the transformation of severely reddened stars, 
the proximity of the reddening path to the data points 
supports the validity of the equations derived by BGB. 
Changing the transformation parameters slightly results 
in a large angle between the data points and the redden- 
ing path. 

A wide spread population of stars is clearly seen along 
the reddening path, as expected from the colour trend dis- 



cussed in Sect. 3.1 (no correction for the varying extinction 
has been applied in this diagram). Again, the stars with 
the lowest reddening within the cluster population are the 
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Fig. 7b is available in jpeg format on astro-ph 



Fig. 7. Ak extinction map, binned in the same manner as the residual map in Fig. || (North is up, East is to the right). 
White spots are positions without stars for evaluation. The individual extinction has been calculated by shifting the 
stars in the K vs. H — K colour-magnitude diagram to a 2 Myr isochrone offset bluewards of the main sequence. 
Transformation to Ak — mag has been performed afterwards, to avoid large errors in the shifting procedure. This 
results in a minimum Ak of 1.86 mag, and a maximum of 4.08 mag, assuming a Rieke & Lebofsky ( |1985[ ) extinction 
law. 
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Fig. 8b is available in jpeg format on astro-ph 



Fig. 8. Ak extinction map derived from HST m205 photometry. See Fig. [?| for details. The coordinate transformed 
HST/NICMOS F205W image is also shown for comparison. Note the different scales (HST/NICMOS: 40" x 40", 
Gemini/Hokupa'a: 20" x 20"). 



bright stars in the Arches cluster center. Moving along the 
reddening line towards higher values of Ay mainly means 
moving radially outwards from the cluster center. As in 
the CMD, a correction for the observed colour gradient 
causes the bulk of the stars to fall onto the main sequence 
with a reddening of Ay — 24 mag, corresponding to the 
cluster center. 



4. Luminosity Functions and incompleteness 
effects 

4.1. Integrated luminosity function 



In Fig. 



11 



we compare the Gemini with the HST lu- 



minosity function for the K^ lans vs. to205 observations. 
For direct comparison of the observational efficiency, no 
colour cut has been applied, but the entire physically rea- 
sonable colour range from approximately to 4 mag in 
H — K', including reddened and foreground objects, has 
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Fig. 10. Two-colour diagram from HST observations. The 
reddening path is shown as a straight line labeled with Ay 
values, and the main sequence is indicated by the thick 
grey line. The area between the dashed lines marks the 
region of reddened main sequence stars. 



been included in the luminosity function (LF). Therefore, 
these LFs are not the ones from which the mass func- 
tions have been derived. The only selection criterion that 
has been applied is a restriction of the photometric un- 
certainty in both magnitude [gk 1 < 0.2 mag) and colour 
{oh-k 1 < 0.28 mag, corresponding to ok' < 0.2 mag 
and an < 0.2 mag). The uncertainty selection in colour 
allowed us to select only those objects that have been de- 
tected with high confidence in both H and K' images in 
the Gemini data, and in the F160W and F205W filters in 
the NICMOS data, respectively. The colour-uncertainty 
selection on the HST data simulates the matching of H 
and K 1 detections used on the Gemini data for the selec- 
tion of real objects. Thus, only objects that are detected in 
both H and K' have been included in the luminosity func- 
tion. This gives us some confidence that we are not picking 
up hot pixels or cosmic ray events. The area covered with 
Gemini has been selected from the HST photometry as 
displayed in Fig. [j] 

As can be seen in Fig. [H], many objects are missed by 
Gemini in the fainter regime, though the actual limiting 
(i.e., cut-off) magnitudes are the same in both datasets. 
This is due to the fact that 50 % of the light is distributed 
into a halo around each star. These halos prevent the de- 
tection of faint objects around bright sources, especially 
in the crowded regions. This effect is most obvious when 
examining the star-subtracted frames resulting from the 
DAOPHOT allstar task. In these frames the cluster cen- 
ter is marked by a diffuse background, enhanced by ~ 20 
counts in K' and ~ 40 counts in H above the observational 
background of 2 and 4 counts in the cluster vicinity, re- 
spectively. In addition to the simple crowding problem due 
to the stellar density affecting both datasets, the overlap 
of many stellar halos hinders the detection of faint stars in 
the Gemini data. At larger radial distances from the clus- 
ter center, more and more faint stars are detected both in 
the HST as well as in the Gemini data (Fig. O). 




14 16 18 

K'trans / m205 



Fig. 11. Comparison of Gemini versus HST luminosity 
functions. 



The fact that the incompleteness corrected Gemini LFs 
follow closely the shape of the HST LFs supports the re- 
sults of our incompleteness calculations, which will be used 
to determine the incompleteness in the mass function. 

4.2. Radial variation of the luminosity function 

Radial luminosity functions were calculated in AR = 5 " 
bins, using the same uncertainty selection as in Fig. 
(Sect. 4.1). The resulting radial LFs are shown in Fig. n2T 
along with the incompleteness determined for each radial 
bin. In the cluster center (lowest panel), the very good 
match of the Gemini and HST LFs for ^rans < 16 mag 
reveals the comparable spatial resolution obtained in both 
data sets. Despite the strong crowding seen already in 
these bright stars, the Gemini AO data resolve the sources 
in the cluster center nicely. When we move on to fainter 
magnitudes, however, we are limited by the halos of these 
bright stars, as discussed above. The clear decrease be- 
low -f^trans = 16 mag marks the point where stars are lost 
due to the enhanced background. When we move radially 
outwards, the limiting magnitude above which faint stars 
are lost shifts towards fainter magnitudes. The tendency 
to loose the faint tail of the magnitude distribution never- 
theless remains clearly seen, though it becomes much less 
pronounced for R > 10", where the Gemini and HST LFs 
resemble each other. For R > 15" (upper panel) we are 
limited by small number statistics due to the small area 
in this radial bin. As it is hard to observe a well-defined 
LF at these radii, we will add the two upper bins when we 
create the radially dependent mass functions in Sect. |5.3[ 
The magnitude-dependent distribution of stars within 
the cluster is evident in these LFs. While bright stars are 
predominantly found in the cluster center, their number 
density strongly decreases with increasing radius. When 
we analyse the Gemini LFs more quantitatively, we find 
25 (50) stars with # t ' rans < 13(14) mag within R < 5", 
but only 8 (23) such stars with 5 < R < 10", and beyond 
10", we observe only 7 (11) such stars. The numbers for 
HST are comparable in the bright magnitude bins. On 
the other hand, the number of faint stars with i^ t ' rans > 
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Fig. 12. Radial variation of the luminosity function. The 
comparison of Gemini and HST LFs is shown together 
with the Gemini incompleteness calculation. The depen- 
dence of the magnitude limits on the distance to the clus- 
ter center is striking. 

18(19) mag increases from 1 (0) to 14 (3) to 48 (16). As 
we see significantly more faint stars in the HST data, the 
corresponding numbers are higher, i.e. the number of stars 
with to205 > 19 mag is in the innermost bin, 24 in the 
intermediate bin, and 81 in the outermost bin. Despite the 
fact that the area on the Gemini frame increases by about 
a factor of 3 between the inner and intermediate bin, the 
number of bright stars is strongly diminished beyond a few 
arcseconds, while the number of the faint stars increases 
by much more than the change in area can account for. 
Although for the fainter stars the effects of crowding and a 
real increase in the fainter population of the cluster cannot 
be disentangled, the decrease in the number of bright stars 
is a clear indication of mass segregation within the Arches 
cluster. 

5. Mass Function 

The mass function (MF) may be defined as the number of 
stars observed in a certain mass bin. The mass function 
in stellar populations is most frequently fitted by a power 
law, whose slope depends on the mass range analysed (e.g, 
Kroupa 2001). In the logarithmic representation, the mass 



d(logiV)/d(logM) - M r 

where M is usually given in solar masses, and T is the slope 
of the mass function. The exact shape of the mass function 
and, in particular, its slope , has been su bject of int ense 
discu ssion (see, e.g., Massey 1995a , 1995b; Scal o 1986|; and 
Scalo |1998| for a review). Massey et al. ( |l995a| , |l995b[ ) re- 
port slopes between —0.7 and —1.7 with a weighted mean 
of —1.1 ±0.1 for M > 7 Mq in several starforming regions 
in the Milky Way, and —1.1 < V < —1.6 with a mean 
of -1.3 ± 0.1 for M > 2b M Q in the LMC. From these 
studies in young clusters and associations, the mass func- 
tion has been suspected to be rather universal, following 
a Salpeter (1955) power law with a slope of V = —1.35 
for stars with masses M > 1 Mq (Kroupa 2001). For the 
two starburst clusters studied in the Milky Way, Arches 



(FKM) and NGC3603 (Eisenhauer |199J), a flat slope of 
~ —0.7 is observed for masses M > 10M Q in Arches, and 
M > IMq in NGC3603. These values have to be kept in 
mind as a comparison for the results presented in the fol- 
lowing sections. The mass regime which will be discussed 
here spans a mass range of 6 < M < 65 Mq for Arches 
cluster stars. 



5.1. Integrated mass function 

The present-day mass function (Fig. |l3|) of the Arches 
cluster has been derived from the colour-magnitude dia- 
gram by transforming stellar luminosities into masses via a 
2 Myr isochronc from the Geneva basic set of stellar evolu- 



function is defined as 



tion models (Lejeune & Schaerer 2001) using the method 
described in Grebel & Chu ( 2000p . Enhanced mass loss 
models were also used, but did not alter the resultant 
mass function. Stellar evolution (mass loss, giant evolu- 
tion) is not important on timescales of the Arches age 
of ~ 2 Myr for stars with initial masses of M < 50 Mq 
(\og(M /Mq) < 1.7, i.e. stellar evolution affects the two 
upper mass bins of the MF at most). No attempt has been 
made to reconstruct the initial mass function (IMF) from 
the present-day MF for M > 50 Mq. A distance modulus 
of 14.5 mag and an extinction of Ay = 24.1 mag have 
been applied. 

The slopes of all mass functions discussed have been 
derived by performing a weighted least-squares fit to the 
number of stars per mass bin. The size of the mass bins was 
chosen to be 5\og(M /Mq) = 0.1 as the best compromise 
between mass function resolution and statistical relevance. 
This bin size is significantly larger than the photometric 
uncertainty in the considered mass and thus magnitude 
ranges. Only mass bins with a completeness factor of > 
75% have been included in the fit. 

Note that we have not attempted to subtract the field 
star contribution. As can be seen in Fig. |l|, the Gemini field 
is mostly restricted to the densest cluster region. When 
comparing to an arbitrary part of the GC field, we do not 
expect to observe the same distribution of stars as in the 
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Arches field, as the faint, reddened background sources are 
negligible due to the high density of bright sources in the 
cluster area. 

In addition, the stellar density in the GC is strongly 
variable, imposing additional uncertainties on the field 
contribution. Neither the Gemini nor the HST field cov- 
ers enough area to estimate the field star population in 
the immediate vicinity of the cluster. A main sequence 
colour cut (1.15 < H — K < 1.90 mag) has been applied 
to the colour-corrected CMDs to select Arches members, 
excluding blue foreground and red background sources. 

To allow for a direct comparison with the results ob- 
tained in FKM, we have used isochrones calculated for a 
metallicity of Z = 0.04 for all MF derivations. The derived 
MFs are displayed in Fig. 13, The overall mass function 
derived from the Gemini data displays the same slope as 
derived from NICMOS within the uncertainties, namely 
r G cmi„i = -0.77 ± 0.16 and r HS T = -0.82 ± 0.14 fitted 
for 10<M<65 Mq (Fig. |r|), which may be extrapolated 
down to 6 Mq when taking into account the incomplete- 
ness correction. The present-day upper mass of 65 Mq 
corresponds to an initial mass of about 100 M© accord- 
ing to the Geneva models. FKM derive an overall slope of 
T = -0.7±0.1 in the inital mass range 6 < M < 120 Mq, in 
good agreement with our results. The remaining difference 
in the maximum mass is due to the different extinction 
and the extinction corrections applied, which represent 
the largest uncertainties in the mass function derivation. 
As in particular the correction of the K magnitude for the 
varying extinction is uncertain due to the unknown extinc- 
tion law of the NICMOS niters, we have also derived the 
mass function for uncorrected K magnitudes, with only 
the colour correction applied, which is independent of the 
extinction law assumed. In this case, the MF appears flat- 
ter with a slope of V ~ —0.5 ± 0.2 (Fig. [r| lower panel). 
The discrepancy in the derived slopes clearly shows that 
the effects of differential extinction arc not negligible, es- 
pecially when deriving mass functions for very young re- 
gions, where the reddening varies significantly. 

Furthermore, we have checked the effect of the binning 
on the MF by shifting the starting point of each bin by one 
tenth of the bin-width, <51og(M/M Q ) = 0.01. The resul- 
tant slopes range from — 0.69± 0.13 < Tccmim < — 0.90± 
0.15 and -0.79±0.12 < r HS T < -0.98±0.13. The average 
slopes for Gemini and HST, Tccmim = — 0.77± 0.15 and 
Thst = —0.86 ± 0.13, respectively, agree well within the 
errors. The slightly flatter slope observed in the Gemini 
data may reflect the more severe incompleteness due to 
crowding. Although all slopes are consistent within the 
errors, the range in slopes derived by scanning the bin- 
step shows that statistical effects due to the binning may 
not be entirely neglected in the MF derivation. 

The metallicity within the immediate Galactic Center 
region has been a matter of discussion during the past 
decade. Several authors report supersolar metallicities de- 
rived from CO index strength and TiO bands in bulge 
sta rs (Fr o gel k Whitford |1987|; Ric h |1988| ; Terndrup et 
al. |l99"0] |i99ll). Carr et al. (2000|) measure [Fe/H] = 



-0.02 ± 0.13 dex for the 7 Myr old supergiant IRS 7, 
and Ramirez et al. (2000) derive an average of [Fe/H] 
= +0.12 ± 0.22 dex for 10 young to intermediate age su- 
pergiants, both very close to the solar value. 

Using a 2 Myr isochrone with solar metallicity Z=0.02, 
the average slope of all bin steps is Tccmim = — 0.84±0.13 
and Thst = —0.91 ± 0.12. The mass function is thus not 
significantly altered when using solar instead of enhanced 
GC metallicity. We note, however, that a lower metallic- 
ity (i.e., in this case solar) steepens the MF slightly, thus 
working into the same direction as the incompleteness cor- 
rection. 

FKM report a flat portion of the MF in the range 
15 < M < 50 Mq, which is not seen in the Gemini MF. This 
plateau can, however, be recovered, when we create a MF 
from _ftT-band magnitudes uncorrected for differential ex- 
tinction, and use a lowest mass of log(Afi ow /MQ) = 0.25. 
For log(Mi ow / Mq) = 0.20 the plateau is seen neither with 
nor without extinction correction. This, again, shows the 
dependence of the shape of the MF on the extinction cor- 
rections applied, as well as on the chosen binning. 

From the considerations above, we conclude that the 
overall mass function of the Arches cluster has a slope of 
T = -0.8 to -0.9 in the range 6<A/<65 Mq. Although 
the uncertainty of missing lower mass stars in the immedi- 
ate cluster center remains, the incompleteness correction 
strongly supports the derived shape of the MF. If the flat 
slope would be solely due to a low recovery rate of low- 
mass stars in the cluster center, this should be visible in 
a much steeper rise of the incompleteness corrected MF 
in contrast to the observed MF. We thus conclude that 
the slope of the MF observed in Arches is flatter than the 
Salpeter slope of T = —1.35, assumed to be a standard 
mass distribution in young star clusters. Such a flat mass 
function is a strong indication of the efficient production 
of high- mass stars in the Arches cluster and the GC envi- 
ronment. 



5.2. Effects of the chosen isochrone, bin size, and 
metallicity 



Blum et al. ( 2001 ) estimate a cluster age of 2 — 4.5 Myr 
for Arches assuming that the observed high-mass stars 
are of type WN7. If we compare the Geneva basic grid 
of isochrones with fundamental parameters obtained for 
WN7 stars (Crowther et al. 1995 ), a reasonable upper age 
limit for this set of isochrones is ~ 3.5 Myr. Crowther com- 
pares the parameters derived for WN stars with evolution- 
ary models at solar metallicity from Schaller et al. (1992) 



and with the mass-luminosity relation for O supergiants 
from Howarth & Prinja (1989), yielding a mass range of 
20 — 55 M©, but with high uncertainties at the low-mass 
end. The more reliable mass estimates for the colour and 
magnitude range observed for WN stars in Arches are re- 
stricted to 35 < M < 55 Mq. From spectroscopic bina- 
ries, the masses of two WN7 stars are determined to be 
- 30 Mq and > 48 Mq (Smith & Maeder |l989t) . The the- 
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oretical lower limit to form Wolf-Rayct stars is 25 Mq for 
the Geneva models (Schaerer et al. 1993). 

In the Geneva basic grid of models with Z = 0.04, the 
3.5 Myr isochrone is limited by a turnoff mass of 32 Mq. 
We have thus calculated mass functions for isochrones 
with ages 2.5, 3.2, and 3.5 Myr in addition to the 2 Myr 
case discussed above. Though the derived slopes scatter 
widely, irrespective of the isochrone used, the slope of the 
MF tends to be even Hatter for any of the older population 
models. We thus conclude that, regardless of the choice of 
model and parameters, the Arches mass function displays 
a flat slope. 

5.3. Radial variation in the mass function 

The radial variation (Fig. |l4|) of the mass function is par- 
ticularly interesting with respect to YC evolution. We 
have analysed the stellar population in Arches within 
three different radial bins, 0" < R < 5", 5" < R < 10" 
and 10" < R < 20". The resulting mass functions for the 
Gemini and HST datasets, along with the radius depen- 
dent incompleteness correction for the Gemini MFs, are 
displayed in Fig. [HJ. 

We confirm the flat mass function observed by FKM 
in the innermost regions of the cluster, which steepens 
rapidly beyond the innermost few arcseconds. Most of the 
bright, massive stars are found in the dense cluster center, 
where the mass function slope is very flat. FKM derive a 
slope of r = —0.1 ± 0.2 from the HST data in this region, 
which is consistent with Fig. O. It is obvious from the 
lowest panel in Fig. [l4] that we are crowding limited in 
the innermost region. We have thus not tried to fit a slope 
for R < 5". 

In the next bin, 5" < R < 10", the mass function 
obtained from the weighted least-squares fit has already 
steepened to a slope of T = —1.0 ± 0.3. Again, the MF in 
this radial bin remains significantly flatter (r = — 0.5±0.4) 
when no A^-correction is applied. Beyond 10" (0.4 pc, 
upper panel) , a power law can only be defined in the range 
10<A/<30 Mq (1.0 < log(M/M Q ) < 1.5), where a slope 
of r = —1.69 ± 0.66 is found, consistent with a Salpeter 
(r = —1.35) law. The large error obviously reflects the 
small number of mass bins used in the fit. Nevertheless, 
Fig. |lj clearly reveals the steepening of the MF very soon 
beyond the cluster center. 

For a more quantitative confirmation of the mass seg- 
regation present in the Arches cluster, we created cumula- 
tive functions for the mass distributions in the three radial 
bins (Fig. |l5|). We have applied a Kolmogorov-Smirnov 
test to quantify the observed differences of these functions. 
When the central, intermediate, and outer radial bin are 
compared pairwise, we obtain in each case a confidence 
level of more than 99 % that the mass distributions do 
not originate from the same distribution. 

Thus, the inner regions of the cluster are indeed skewed 
towards higher masses either by sinking of the high mass 
stars towards the cluster center due to dynamical pro- 



cesses or by primordial mass segregation, or both. The 
same effect is observed in the similarly young cluster 
NGC3603 (Grebel et al. 2002, in prep.). 

5.4. Formation locus of massive stars 

We find several (5) high mass stars in the cluster vicin- 
ity. These stars fall onto the Arches mai n se quence after 
applying the reddening correction (Sect. 3J). When the 
entire HST field is separated into two equalsize areas, the 
first area being a circle of radius 16" around the cluster 
center, and the second the surrounding field, no additional 
comparably high-mass stars are found with Arches main 
sequence colours except for the two bright stars found at 
the edges of the G emini field. In the dynamical models of 
Bonnell & Davies ( |1998| ) there is a low, but non-zero prob- 
ability that a massive star originating outside the clus- 
ter's half-mass radius might remain in the cluster vicinity. 
High-mass stars formed near the center show, however, 
a tendency to migrate closer to the cluster center. The 
disruptive GC potential might enhance the ejection pro- 
cess of low-mass stars, but according to equipartition, it 
is unlikely that the most massive stars gain energy due to 
dynamical interaction with lower mass objects. We have 
to bear in mind, however, that interactions between mas- 
sive stars in the dense cluster core could cause the ejection 
of high-mass stars. 

N-body simulations performed by Portegies Zwart et 
al. (1999) show that the inclusion of dynamical mass segre- 
gation in cluster evolution models enhances the collision 
rate by about a factor of 10 as compared to theoretical 
cross section considerations. For a cluster with 12000 stars 



initially distributed according to a Scalo ( pL986[) mass func- 
tion, a relaxation time of 10 Myr, and a central density 
and half-mass radius comparable to the values found in 
Arches, about 15 merging collisions occur within the first 



10 Myr (lt r , 



of the simulated cluster. Shortly after 



the start of the simulations, frequent binary and multiple 
systems form from dynamical interactions leading to the 
ejection of several contributing massive stars. The flat MF 
in the Arches center as compared to a Scalo MF, contain- 
ing a larger fraction of massive stars to interact, may even 
increase the collision rate. 

Though it is likely that the high-mass stars seen in 
the immediate vicinity of the Arches cluster have formed 
from the same molecular cloud at the same time as the 
cluster, a final conclusion on the possible ejection of these 
stars from the cluster core due to dynamical processes can 
only be drawn when velocities for these cluster member 
candidates are available. 



6. Comparison with cluster formation models 

In this section, we will first use the simple analytical ap- 
proach to analyse internal cluster dynamics as summarised 
in Binney & Trcmaine (1987), and will then compare the 
derived dynamical timescales to N-body simulations where 
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Fig. 13. Arches mass function derived from the 
Gemini/Hokupa'a colour-magnitude diagram shown 
in Fig. ^ A 2 Myr main sequence isochrone from the 
Geneva set of models (Lejeune & Schaerer 2001 ) was used 
to transform magnitudes into stellar masses. The mass 
function has been derived for bins of <51og(M/M ) = 0.1 
with a lowest mass bin log(M/M Q ) = 0.20. 
Upper panel: with Ak correction 
Lower panel: without Ax correction 



the tidal force of the Galactic Center potential is consid- 
ered. 

From the transformation of stellar magnitudes into 
masses, a rough estimate of the timescales relevant for 
dynamical cluster evolution can be made. For larger area 
coverage, we have used the HST/NICMOS data in the 
calculation below. The timescales characterising cluster 
evolution are the cluster's crossing time, t CTOSS , which is 
simply given by some characteristic radius divided by the 
average velocity, i.e., the mean velocity dispersion of the 
cluster, r c /(a), and the relaxation time, i r h- The median 
relaxation time is the time after which gravitational en- 
counters of stars have caused the system to equilibrate to 
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Fig. 14. Radial change of the mass function as observed 
in the Gemini/Hokupa'a data. A very flat mass function 
is seen in the inner cluster regions, where predominantly 
massive stars are found. The slope of the mass function 
increases towards the Salpeter value (T = —1.35) already 
at a radial distance from the cluster center of > 10" (0.4 
pc at a distance of 8 kpc). 

a state independent of the original stellar orbits (Binney 
& Tremaine 



1987, hereafter BT87) 



irh 



M 



6.63 • 10 8 yr 
ln(0.4AT) V !0 5 M G 



1/2 



1 Mr. 



1 pc 



3/2 



where M is the total mass within some characteristic ra- 
dius r c , m* is a characteristic stellar mass, here defined as 
the median of the observed mass distribution, and N is the 
number of stars in the cluster. i r h may also be expressed 
in terms of the crossing time, i r h ~ 0.1 N/ In AT ■ t CIOSS . 

To evaluate the above formula, we need to know the 
total mass of the system, M, a characteristic radius, r c , 
the number of stars, N, and the characteristic stellar mass, 
m*. In principle, we are able to derive most of these quan- 
tities from isochrone fitting, by assigning each star a mass 
corresponding to its if-band luminosity and analysing the 
resulting spatial distribution of masses. Naturally these 
individual masses cannot be accurate for each star, as they 
depend strongly on the choice of the isochrone and inherit 
the uncertainties of the photometry. The integrated prop- 
erties are, however, not very sensitive to the age of the 
model isochrone used, within the reasonable age range 
of Arches, ~ 2 — 3 Myr (results will be given below). 
We have attempted to create a density profile from the 
HST/NICMOS data, but as the profile appears to be very 
distorted, we have chosen to use the half-mass radius, r nm , 
as a characteristic scale, rhm has been obtained from the 
observed stellar population of the cluster under the as- 
sumption that the mass in the cluster center, relevant for 
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Fig. 15. Cumulative functions of the stars in each radial 
bin obtained from HST/NICMOS data. Each radial bin 
corresponds to one radial mass function in Fig. [l4| The 
cumulative distributions have been normalised at 10 M Q , 
the limit where all radii are more than 75 % complete. 



the spatial scale on which gravitational interactions are 
important, is dominated by the detected high-mass stel- 
lar population in the cluster center. The relaxation time 
derived from rh m is usually referred to as the "half-mass 



relaxation time" (cf. Portegies Zwart et al. 2002). 

The mass distribution within the cluster as derived 
from the isochrone allows us to estimate rhm, M, and m*. 
Taking all stars within a certain radius with Arches main 
sequence colours as cluster members, we can also estimate 
N. 

We obtain a half-mass radius of rh m = 10" = 0.4 pc, 
irrespective of the age of the isochrone used (2 Myr or 
3.2 Myr) to transform A"-band magnitudes into masses. 
The total mass measured within this radius is -M(r mn ) = 
6300 M @ (5400 M Q for the 3.2 Myr isochrone). By extra- 
polation of the mass function down to 0.1 M Q , FKM de- 
rived a total mass of 12000 Mq within 9" radius. To esti- 
mate the total mass in the cluster center, they used the 
relative fraction of stars with M > 40 Mq in two annuli 
separated at R = 3" as a scaling factor. This procedure 
ignores the segregation of high-mass stars. As the clus- 
ter shows strong evidence for mass segregation (see Sect. 
5.3), this may overestimate the total cluster mass (within 
R < 9"). On the other hand, our mass estimate, ignor- 
ing incompleteness and the lower mass tail, is a lower 
limit to the total mass within R < 10". We thus con- 
clude that for an order-of-magnitude estimation, a total 
mass of 10 4 Mq within the half-mass radius is a good 
approximation. This yields an average mass density of 
p(.R<r U m)~4-10 4 M Q pc- 3 . 

The calculation of the relaxation time requires an es- 
timate of the characteristic stellar mass and the number 
of stars. The median mass within R < 10" is 7.3 M Q in 
our observed mass range of 1.4 < M < 65 Mq in the HST 



sample, and the number of stars actually observed is 486. 
We choose to use the median mass here as a more realistic 
mass estimate than the mean for individual stars partic- 
ipating in interaction processes, as we are aware of the 
fact that the mean is highly biased by the high fraction 
of high-mass stars in the cluster center. Nevertheless, we 
have to keep in mind that this median value is still unre- 
alistically high as we miss the low-mass tail of the mass 
distribution, which also has the highest number of stars. 

As the low-mass tail of the distribution is highly incom- 
plete, we estimate N to be at least 10000 M©/7.3 Mq ~ 
1400. This results in = 1.2 Myr. Decreasing the charac- 
teristic mass and increasing the number of stars lengthens 
the relaxation time. For instance, lowering m* to 1-Mq, 
where we take into account the fact that the cluster cen- 
ter is mass segregated as opposed to a standard Salpeter 
mass function and thus the characteristic mass of a star 
should be higher, and raising N to 10000 stars, would re- 
sult in £ r h = 6.2 Myr. If a significant fraction of low-mass 
stars would have been formed and survived in the cluster 
center, the half-mass relaxation time would increase even 
more. The present-day relaxation time is thus at least on 
the order of or longer than the age of the Arches cluster. 

A relaxation time longer than the lifetime of the cluster 
would imply that not yet sufficient time has passed for dy- 
namical mass segregation to take place and would thus in- 
dicate that primordial mass segregation was present at the 
time of cluster formation. However, we have so far ignored 
the Galactic Center tidal forces accelerating the dynami- 
cal evolution of the cluster. Kim et al. (2000) performed 
N-body simulations using the observed parameters of the 
Arches cluster to trace the dynamical evolution of the clus- 
ter and constrain initial conditions. They use a total mass 
of 2 • 10 4 M Q , tidal radius of ~ 1 pc and a single power-law 
IMF with slopes of T = -0.5, -0.75,-1.0 and -1.35. The 
number of stars ranges between 2600 for a lower mass cut- 
off of m low = 0.1 Mq and 12700 for m\ ovl = 1 Mq. From 
comparison with the HST/NICMOS data of FKM their 
models yield a power law with T = —0.75 as the most 
probable initial mass distribution. Although this slope is 
very close to the observed present-day MF, they also note 
that mass segregation takes place on timescales as short as 
1 Myr, such that the cluster looses all memory of the initial 
conditions shortly after formation. Portegies Zwart et al. 
(2002) use a Scalo (1986) MF to model the Arches cluster, 
and derive a total cluster mass of 4- 10 4 M Q , and an initial 
relaxation time of 20-40 Myr. When comparing with the 
same set of HST/NICMOS observations, they conclude 
that a standard IMF can evolve into the current MF due 
to dynamical segregation, and that the IMF did not need 
to be overpopulated in high-mass stars. Their computa- 
tions indicate that the relaxation time is strongly variable. 
During the dynamical evolution, t r h increases strongly af- 
ter core collapse, which occurs within ~ 2 Myr, due to 
cluster re-expansion, and only starts to decrease after 10 
Myr or later. The observed relaxation time does thus prob- 
ably not trace the cluster's initial conditions, but reflects 
the current dynamical state in the evolution of Arches. In 
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particular, a present-day relaxation time larger than the 
cluster's age does not necessarily imply that the cluster is 
not dynamically relaxed. 

Unfortunately, this means that we are not able to dis- 
tinguish between primordial and dynamical mass segrega- 
tion from the estimated relaxation time. For the Orion 
Nuclear Cluster (ONC) and its core, the Trapezium, 
which has been studied in greater detail (e.g., Hillenbrand 
1997, Hillenbrand k Hartmann 1998), Bonnell & Davies 
(1998) derive dynamical timescales too long to explain 
the segregation observed in high-mass stars within the 
Trapezium by dynamical evolution, concluding that a sig- 
nificant amount of primordial segregation must have been 
present. In the case of the Arches cluster, the external 
gravitational field acts towards a fast disruption, thereby 
impeding the equilibration process, such that dynamical 
segregation may well be under way. 

As the different model calculations do not agree 
with respect to the IMF required to create the ob- 
served present-day mass distribution, a final conclusion 
on whether or not the Arches initial mass function had to 
be enriched in massive stars, thus supporting high-mass 
star formation models, cannot be drawn. 

The evaporation time, setting the scale for dynami- 
cal dissolution by internal processes, can be estimated as 
icvap » 136 Wax ~ 177 Myr (857 Myr for m* = 1M , 
BT87). Kim et al. ( |1999| ) have shown that the time re- 
quired to disrupt an Arches-like cluster within the GC 
potential is only 10 Myr, much shorter than evaporation 
by internal processes will ever be relevant. The models 
by Portegies Zwart et al. ( |2002 ) suggest somewhat longer 
evaporation times in the range 30 < t cva p < 50 Myr for an 
Arches-like cluster. The external potential is thus the dom- 
inating factor in the dynamical evolution of Arches. Note 
that, however, even a cluster as dense as Arches would 
not survive for more than 1 Gyr independent of its lo- 
cus of formation. The relaxation and evaporation times 
found for Arches are very similar to our results for the 
young, compact cluster in NGC3603 (Grebel et al. 2002, 
in prep.). However, as NGC3603 is not torn apart by ad- 
ditional external tidal forces, it may survive much longer 
than Arches. For comparison, massive young clusters in 
the Magellanic Clouds have relaxation times of ~ 10 8 yr, 
and corresponding evaporation times of ~ 10 10 yr (e.g. 
Subramaniam et al. 



1993), and may thus survive for one 



Hubble time. 



7. Conclusions 

We have analysed high- resolution Gcmini/Hokupa'a adap- 
tive optics and HST/NICMOS data of the Arches cluster 
near the Galactic Center with respect to spatial varia- 
tions in the mass function and their implications for clus- 
ter formation. A detailed comparison of the Gemini data 
to HST/NICMOS observations allows us to investigate 
the instrumental characteristics of PSF fitting photom- 
etry with the Hokupa'a AO system. 



7.1. Technical comparison of Gemini '/Hokupa'a with 
HST/NICMOS 

The calibration of the Gemini/Hokupa'a data of the 
Arc hes cl uster using HST/NICMOS data from Figer et 
al. ( 199E ) allows us to carry out a detailed technical 
comparison of the two datasets. Maps of photometric 
residuals show a strong dependence of the calibration er- 
ror on the stellar density within the field. In particular, 
the vicinity of fainter objects to bright stars causes the 
Gemini magnitude to be underestimated in comparison 
with HST/NICMOS. This is understandable as the un- 
compensated seeing halos of bright stars enhance the back- 
ground in a non-homogeneous way, thereby causing an 
overestimation of the background and a subsequent under- 
estimation of the faint objects' magnitude. Conversely, the 
flux of very bright sources seems to be overestimated. The 
correlation of the photometric residual with the position of 
bright stars is more pronounced in K' , where crowding is 
the dominant source of photometric uncertainty, while the 
effect of angular anisoplanatism is less severe. In the H- 
band the anisoplanatism is more pronounced, and hence 
photometric uncertainties are a blend of uncertainties due 
to the distance to the guide star and due to the proximity 
to bright sources. 

The incompleteness of the luminosity function mea- 
sured with Hokupa'a increases faster at fainter magnitudes 
than the incompleteness observed in the NICMOS LF 
despite the comparable detection limit in both datasets. 
As expected, this effect is particularly pronounced in the 
dense cluster center, where crowding is most severe. 

As the Strehl ratio determines the amount of light scat- 
tered into the seeing induced halo around each star, a good 
SR is crucial to achieve not only a diffraction limited spa- 
tial resolution, but to benefit from the adaptive optics 
correction in dense fields containing a wide range of mag- 
nitudes. In the Arches dataset, the SR of only 2.5 % in 
H and 7 % in K' as compared to 95 % in F160W and 
90 % in F205W (NICMOS) is clearly the limiting factor 
for crowded field photometry. As Hokupa'a was initially 
designed and developed for the 3.6m CFH telescope, its 
performance at the 8m Gemini telescope is naturally con- 
strained by the limited number of only 36 actuators. In the 
case of the Arches science demonstration data, additional 
constraints were given by the seeing, the high airmass due 
to the low latitude of the Galactic Center, and the guide 
star magnitude. Under better observing conditions and 
with a brighter guide star Strehl ratios of up to 30 % can 
be achieved with Hokupa'a at Gemini. Higher order AO 
systems are currently capable to produce SRs of up to 
50%. 

The Gemini ground-based AO data are comparable to 
the HST/NICMOS data in the resolution of bright sources 
(K' < 18 mag), and in a non-crowded field. They do reach 
their limitations in the densest cluster area and in the 
case where faint stars are located close to a bright object. 
Higher Strehl ratios would of course reduce this unequal- 
ly- 
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7.2. Photometric results 



7.4. Cluster dynamics 



A strong colour gradient is detected over the field of the 
Arches cluster, revealing an increase in visual extinction 
of approximately 9 < A Ay < 15 mag when progressing 
outwards from the cluster center. The visual extinction is 



estimated from the Rieke & Lebofsky ( 1985 ) extinction 
law to be Ay ~ 24 mag in the cluster center, increasing 
to a maximum of 33 < Ay < 39 mag in the vicinity, in ac- 
cordance with Cotera et al. (2000), who found a maximum 
of Ay — 37 mag in the Arches field. Within the central 
5" radius, however, no colour gradient is observed. This 
indicates that the cluster center has been stripped of the 
remaining dust and gas either by strong stellar winds from 
massive stars or by photo-evaporation, or both. Beyond 
5", a linear increase in extinction is observed, suggesting 
an increasing amount of dust with distance from the clus- 
ter center. Photo-evaporation due to the strong UV ra- 
diation of the 8 WN7/8-stars and more than 100 O-stars 
found in the cluster center is most probably responsible 
for dust dissolution. 

The TO160 - to205, to205 (equivalent to H - K, K) 
colour-magnitude diagrams derived from the HST and 
Gemini data sets both show a bent main sequence fol- 
lowing this colour trend. The main sequence straightens 
out when correcting for this colour variation. A spatial 
analysis of the CMDs reveals the bulk of the bright stars 
on the Arches field to be located in the cluster center. 



7.3. Mass functions 

Present-day mass functions have been derived from the 
CMDs after linear correction of the colour trend and corre- 
sponding change in extinction over the field, and selection 
of a reasonable main sequence colour cut. The integrated 
mass function derived from the Gemini photometry dis- 
plays a slope of r = -0.8 ± 0.15 for 6 <M< 65 M , less 
steep than the Salpeter slope of T = —1.35. This value 
agrees with the slopes derived from the HST data in the 
same manner (Sect. |5|), and with the values presented in 
FKM within the uncertainties. When the magnitudes are 
not corrected for differential extinction, the slope of the 
MF is significantly flatter, T ~ —0.5 ± 0.2. Particularly 
in young star forming regions, the effects of differential 
extinction are thus clearly non-negligible. 

The analysis of the radial dependence of the mass func- 
tion reveals a very flat IMF in the immediate cluster cen- 
ter with a slope close to zero. The IMF slope seems to 
increase outwards with T = —1.0 ± 0.3 for 5 < R < 10" 
and T = -1.7 ± 0.7 beyond R > 10". We have created 
cumulative functions for the stars in each radial bin, and 
performed a KS test to derive the significance level of the 
variance in the mass distributions. The probability for the 
observed distributions to originate in the same mass func- 
tion is below 1% when comparing each two of the three 
radial bins analysed. 



The flat mass function in the Arches center is a strong 
indication for mass segregation. While the center seems 
to be dominated by high-mass stars, the cluster edges dis- 
play the standard behaviour of a young stellar population. 
A similar radial dependence of the mass function is ob- 
served in the young, compact cluster NGC3603 (Grebel 
et al. 2002, in prep.), located in a normal star forming en- 
vironment in the Carina spiral arm. The fact that two out 
of three compact young clusters found in the Milky Way, 
which have been analysed in such detail to date, display 
a flat mass function slope in the core indicates that such 
a behaviour might be typical for starburst clusters and is 
not restricted to the extreme GC environment. 

Mass segregation in a compact, massive cluster can 
either be caused by an enhanced production efficiency of 
massive stars, or by dynamical segregation during the clus- 
ter's evolution. We roughly estimate the present-day re- 
laxation time of Arches from mass function considerations 
to be about a few Myr (~ 20t cross ), and thus of the same 
order of magnitude as the cluster age of ~ 2 Myr. Again, 
a similar timescale has been found for NGC 3603 as well 
(Grebel et al. 2002, in prep.). The comparison with N- 
body simulations suggests, however, that the dynamical 
evolution of massive clusters close to the Galactic Center 
whipes out the initial conditions within less than 1 Myr. 
We are therefore not able to distinguish between primor- 
dial and dynamical mass segregation. Both effects are in- 
tertwined at the current state of cluster evolution, such 
that a detailed dynamical analysis is crucial for a thorough 
understanding of the formation process. This analysis has 
to await deep, high resolution infrared spectroscopy to ob- 
tain radial velocities and proper motions for a significant 
fraction of stars belonging to the cluster population. 
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